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Abstract 

Using a simple relativistic QFT model of scalar fields we demon- 
strate that the analytic confinement (propagator is an entire function 
[^ I in the complex p^-plane) and the weak coupling constant lead to the 

CN ■ Regge behaviour of the two-particle bound states. 

In QCD we assume that the gluon vacuum is realized by the self- 
^^ I dual homogeneous classical field which is the solution of the Yang- 

^P ' Mills equations. This assumption leads to analytical confinement of 

quarks and gluons. We extract the colorless 0^^ two- gluon state 
from the QCD generating functional in the one-gluon exchange ap- 
proximation. The mass of this bound state is defined by the Bethe- 
Salpeter equation. The glueball mass varies in the region 1470Mev < 
r~| , Mq++ < 1600Mei' for 0.2 < a^ < 0.5 if the gluon condensate is 

^(G^.Gp= 0.012(^6^)4. 
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1 Introduction 

In this paper we shall calculate the mass of the scalar glueball O^^ . Before 
entering into a discussion, formulation and solution of this problem we would 
like to make some general remarks. 

At the present time, QCD is considered to be an uniquely true theory of 
strong interactions, but there are different directions of its realization: lat- 
tice calculations, QCD sum rules, Nambu-Yona-Lasinio model plus or minus 
instantons, chiral perturbation theory, calculation based on the Wilson loop, 
bag and quasipotential models, models based on the analytic confinement, 
and maybe something else. It is useless to give references here because the 
literature is so numerous that it is practicaly impossible to list all the papers 
without missing some of them. 

In our opinion, theoretical reasons for this variety consist in the following. 
From a general point of view the QFT methods can be used if the propagators 
of particles and vertices, which describe interactions, are known and coupling 
constants are weak. Then perturbation calculations are valid and the Bethe- 
Salpeter equation is an adequate relativistic mathematical tool for computing 
bound states in QFT. The quantum electrodynamics, where the propagators 
and vertices are given and the coupling constant a = -^ is small, is the best 
example of this type of calculations. 

However, in the case of QCD the situation is more difficult. Although the 
QCD coupling constant is considered to be relatively small a^ < 1, we can 
not be convinced that propagators of quarks and gluons have the standard 
Dirac and Klein-Gordon form for large distances, i.e., in the hadronization 
region. Unfortunately, a rigorous generally accepted analytic solution of this 
problem in QCD is still missing. Nevertheless, our hope is that a simple 
formalism similar to QED should exist in QCD, but a background form of 
these "basic bricks" is not yet found in QCD. 

Another problem is that we do not have a self-consistent conventional 
theory of bound states in QFT at all and in QCD in particular, which can be 
applied in any cases like the Schrodinger equation is applied to any nonrela- 
tivistic quantum phenomena. Therefore different approaches with a doubtful 
basis are used to calculate bound states in relativistic cases. One can see that 
in most cases authors of different approaches strive to reduce the QFT bound 
state problem to a nonrelativistic Schrodinger equation with an appropriate 
potential. 



The generally accepted physical picture is that processes of hadronization 
take place at large distances where the confinement of quarks and gluons 
plays the main role. It should lead to the following statements. First, the free 
Dirac and Klein-Gordon equations for free particles are not adequate tools 
to describe quarks and gluons in this region. Second, currents, or vertices, 
which describe the connection of quarks and gluons with hadrons, cannot 
be local. In other words, only the knowledge of the vacuum QCD at these 
distances should dictate what equations have to describe quarks and gluons 
and nonlocal vertices. In addition, the QCD effective coupling constant as 
should be small enough in this hadronization region if we want to use the 
Bethe-Salpeter type equations. 

We think that a successful guess about the structure of QCD vacuum at 
large distances can radically solve problems of confinement and hadroniza- 
tion. 

In this paper, we want to show, using a simple model of scalar field, that 
in the frame work of QFT that analytic confinement of constituent particles 
(propagators are entire analytical functions in the momentum p^-space) and 
weak coupling constant (the Bethe-Salpeter equation can be used) lead to 
the Regge spectrum of bound states. 

Thus, if the QCD vacuum results in analytic confinement of quarks and 
gluons and the QCD coupling constant as is small, hadron bound states are 
expected to have the Regge type spectrum and we are able to calculate the 
masses of these states using the Bethe-Salpeter equation. 

In this paper, we calculate the mass of a glueball with quantum numbers 
O"*"^. Glueballs are very interesting objects from both theoretical and exper- 
imental points of view. The experimental situation is not clear up to now 
(see, for example, ||l|, H), and we shall deal with it. However, the theoreti- 
cal problem is quite intriguing: How can massless gluons be joined together 
into massive objects owing to the gluonic QCD self-interaction? It is evi- 
dent that the solution of the glueball problem is directly connected with our 
understanding of nonperturbative Yang-Mills dynamics. Theoretical interest 
in glueballs arised just from the beginning of QCD and undoubtedly will 
continue until the structure of QCD vacuum is understood. All theoretical 
approaches (see, for example, the last papers 0, Q where the previous publi- 
cations are listed) are based on the understanding that the structure of QCD 
vacuum plays the decisive role in the solution of this problem, but all authors 
avoid to make any direct guesses about its explicit structure and take it into 



account in an indirect way only. 

Our approach is based on the assumption that the QCD vacuum is real- 
ized by the self-dual homogeneous vacuum gluon field which is the classical 
solution of the Yang- Mills equations and in the one-loop approximation the 
QCD vacuum energy has minimum for nontrivial configurations of this field. 
This field leads to the analytic confinement of quarks and gluons, so that we 
could expect reasonable values for the glueball mass. This vacuum gluon field 
has been applied to quarkonia in papers P, H, 0]- We shall use the functional 
integral techniques to calculate mass spectrum of bound states (see [^). In 
this techniques the Bethe-Salpeter equation is written in a symmetrized form 
permitting variational calculations. 

2 The toy model 

In this section, we consider a simple quantum field model with confinement 
and demonstrate the features of possible bound states which can be inter- 
preted as standard physical particles. We consider the scalar field ^{x) with 
the Lagrangian in the Euclidean metrics 

C{x) = —-Lp{x)'m^e'^ Lp{x) — gLp^(x). (1) 

The equation for the free field 0(x) looks like 



m 



e'^ip{x) =0, or me '^{p{p) = (2) 



and the solution (f){x) = because the function e~^ 7^ 0, i.e., it has no 
zeroes for any real or complex p"^. Exactly this property means analytic 
confinement. The field (f){x) can exist in the virtual state only, so that it can 
be called virion field (see [H). In addition, one can say that the field (f){x) 
describes constituent particles. 
The propagator 



D{p') = —, e-i. Dix) = fj^^ D{p) e^^- = ^ 
m^ J [271)* [4tt)- 



(3) 



is an entire analytical function in the complex p^-plane. It guarantees the 
confinement of "particles" <^(x) in each perturbation order over the dimen- 
sionless coupling constant 

KAnmJ ~ 
The scale of the confinement region can be characterized by the value 

2 1 Jdx x'^D{x) 8 



^"^ A^ . jdx D{x) 



m 



The mechanism of arising of bound states can be described in the follow- 
ing way. Let us consider the partition function 



5ip e '''''= 5ip e 
The effective Lagrangian with the one-particle exchange reads 

^eff = -li^D-'v) + ^^{^'D^') + .... (5) 

The term {ip^Dip^) can be transformed as 
L2 = // dxidx2 ip'^{xi)D{xi - X2)(P>'^{X2) 

= ^— // dxidx2 (p{xi)(p{x2)D{xi - X2)(p{x2)(p{xi) 

(M! 

2 
where Xi = x + ^y, X2 = x — ^y and 

J(x, y) = (p(x+ -y] (filx- -y] = (p{x) e^ ^v'(a;), 

J+{x,y) = J{x,-y). 

Let us introduce an orthonormal system {UQ{y)}, where Q = {n,l, {/i}) 
are four-dimensional quantum numbers so that 

[dyUQ{y)UQiy) = 5QQ,, and Y.UQ{y)UQiy') = 5iy - y'). (6) 



dx dyi^D{yi) dy2^ D{y2)J{x,yi)5{yi - y2)J^{x,y2), 



Then 



L2 = -^ jdx Jq{x) ■ Jq{x), 
Jq{x) = i^{x)VQ(dMx), J+{x) = '^{x)Vq{- 1))^{x), 
VQid) = fdy./Di^)UQ{y)e^^, 



and 



e ^2 (v^i:'^'') = gxp < ^^4^ > I dx JA[x)Jn[x 



daBB+ exp <^ — = ^yrfx (5+(a;)jQ(a;) + Bq{x)J+{x)) 

= jdaBB^ exp \^Y.jd^ <t^{x)[Bi{x)VQ(d) + i?Q(a;)\/Q(- 9)]0(a;) 
where 

Q 
We substitute this representation into the partition function Z and after 
integration over ip we get 

.. -Y.{B+BQ)~\^r\n{l+?>gV2(B+VQ+BQV+)D+...) 

Z = l[5BQ5B^e ^ 
•^ Q 

- E {Bi[^QQ'"^^QQ']BQ)+0{B^) 

\{5BQ5B^e ««' . (7) 

Q 

Here the polarization operator is 

mQQ^z) = (3(?)^Tr {V^DVq,D) 

dyidy2 UQ{yi)\U{z;yi,y2)UQ'{y2), 



yi -y2\ r^f vi -y2 



XU{z; y,,y2) = {?>gYJD{y,)D z + '-^-^ ]D{z- '-^-^ JD^y^). 



where z = Xi — X2- 

The Fourier transform looks hke 



AnQQ'(p) = Xjdz e'P'UQQiz) = jjdyidy2 UQ{yi)Sp{yi,y2)UQ,{y2), 
whith 

5p(l/i,Z/2) = ^^e-i& ■ e-^(^i-^^^^+^i). (8) 

The orthonormal system {Uq{x)} should diagonalize the kernel Sp{yi,y2) 

In other words, we have to find the spectrum and eigen functions of the 
problem 

SpUq = \q{p^)Uq or J dy2 Sp{yi,y2)UQ{y2) = \Q{p'^)UQ{yi). (9) 

It should be stressed that this equation is nothing else but the Bethe-Salpeter 
equation in the one-boson exchange. In order to go to the standard form of 
the Bethe-Salpeter equation, we have to introduce new function Uq^y) = 

J D{y)'^Q{y) and go to the momentum representation (see, for example, 

i)- 

If we are able to do it, the partition function Z in (0) looks like 

Z= mSBgSB+e « (10) 

Q 

and the equation 

1-Aq(-M2) = (11) 

defines the mass of the bound states with quantum numbers Q. Then 
1 - Aq(p') = -A'q(-M5)(/ + Mj) + Oiip' + Mlf) 

and if 

-A' (-M2) > 



we can interpret the field Bq as a physical particle with quantum numbers 
Q and mass Mq. 

The Bethe-Salpeter equation (^) can be solved and the eigen functions 
look like 



UQ{y) = U^i{,}{y) = CniTi{^}{ay)L'+\aV)e--'' ' 



2 2 

a = m 



Vs 



Here, the angle polynomials Tis^^^{y) = Ti^^^^ ^^^y{y) are normalized by 

Y,Ti{f,}{x)Ti{^}{y) = Cliixy)) 

where C/(t) are the Gegenbauer polynomials. The eigenvalues are 

^ -A 



Ki{p 



{2 + V3) 



l+l+2n 



(12) 



The equation { ^A\ ) gives the spectrum 

Ml = 2mMn [ y J + (Z + 2n) ■ 2m^ ln(2 + ^3) (13) 



nl 



where Ac = (2 + -\/3)^ = 13.928... . One can see that this spectrum has the 
pure Regge behaviour. 
Now we can represent 



i~\q{p' 



X 



^Q 

ei^ 



e 2m^ 



Xc (2 + v^)'+2" 

ZQ(/ + Mj) + 0((p2 + Mj)2) 



with 



A 






1 



^ Ac (2 + v^)'+2- 2m2' 
Thus, the kinetic term in Z (ITO|) is 



{b+{p)[1-Xq{p')]Bq{p) 

= Zq (bUv) [(/ + Ml) + 0((p2 + M^f)] BM) . 



and after the renorinalization 

becomes of the standard form. 

In particular, the lowest scalar bound state Bq has the mass 



^(00) = "^y^lny, 

For example, for A = .1 we get 

M(oo) = 3.142... -m. 
Finally, we conclude: 

• Bound states exist for small coupling constants 

A < Ac ^ 13 
and their masses grow when A — > as 



Mq ~ mW2 In y. 

Analytic confinement leads to the pure Regge spectrum for all bound 
states with quantum numbers Q = (nl) and the slope of the Regge 
trajectories is defined by the scale of the confinement region A^„j and 
does not depend on the coupling constant A. 

We can see that the size of the confinement region Vcnf @ excels re- 
markably the Compton length of all bound states 

It means that the physical particles which are described by the fields 
Bq{x) and all transformations with them take place in the confinement 
region induced by the constituent particles (f){x). 



If A < .1, i.e., we have the standard condition for perturbation calcu- 
lations, then 



M > m^2 In(lOAc) = 3.142... ■ m. 

In other words, M remarkably excels the "initial mass" m. If A is small 
enough, the bound state can be quite large. 

This scheme will be applied in main features to find the glueball mass in 
QCD (see Table I). 
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Toy model 


QCD 


Gluon propagator is 


QCD vacuum is realized 


assumed to be an entire 


by self dual homogeneous 


analytic function 


gluon field B^ 


e~^' 


Gluon propagator is 




an entire analytic function 


i 
Li = g(j)^{x) 

i 


i 
Lj = gTT{F^,[A^, A,]) 

i 


i 


S^TT{[A^A^]^^Df,^^a(s{xi,X2)[AaA^]^,^) 

I 


^J{Xi,X2)D{xi - X2) J(X2, Xi), 


Firz transformation^ scalar glueball 




^J{Xl,X2))A^{xi - X2)J{X2, Xi) 


J{x,y)=<P{x + y/2)<p{x-y/2) 
= EJQ{x)UQ{y)=EJQ{x)UQ{y) 

Q Q 


Jix, y) = {A^{x + y/2)e'^^-^yU,{x - y/2)) 


Jq(x)=0(x)Vq(H)0(x) 

\ 


Jq{x) = {A,{x)Vq{V)A,{x)) 


Li = 


^EjdxJ^Qix) 
Q 

i 

^' and Integration over (f){x) and A^{x) 
I 

i 

alpeter equation: 

= {UqUUq,) = \q{p^)6qq, 

i 
rum: 1 = Aq(-M^) 


g^(4)^j^5^g-(B|)+s{BQJc 


JDBq e-(-SQ)-Trln[l+9{BQVQ)D] _ 


Bethe-S 
g^TTiVqDVQ'D) 


Mass spect 
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3 QCD Lagrangian in the self-dual gluon field 

Now let us go to QCD. The action has the form 

S=-^ld'xTiG^,G,,, (14) 

Gfj,i/ = o^Ay — OijA^ — tg[A^, A^\, A^ = t A^. 

We suppose (see P, §, 0]) that the QCD vacuum is reahzed by a homogeneous 
(anti-)self-dual gluon field which has the form 

Bf,{x) = B^yXi, = hBh^{x) = fiBh^yXu , (15) 



n = en", 


n+ 


= h, 


(nn) = n'^n'' = 1, 


"^p"pu 


~^HV) 


bpu 


= ^fiuapOa/B = iO^i/- 



After substitution 

AJx) -^ AJx) + BJx) 



one has 



Here 



Gf,u = :F^u-ig[A^,M-2B^ 



flU- 



where we have used the notation 

v;:" = vi^ix) = 5'''^ - in"\{x), v^ = l-- tnb,,x,, 

[V^, VJ = 2inbf,^, n = gBn, fi"^ = ^5n"^ 

The QCD action (|T^ in the background field B looks like 

S = -^Jd^xTi {:F^,T^, + 4tgB^, [A^, A,] 
- 2igP^,[A^, A,] + AB^,B^, - g''[A^,, A,f] . 

12 



In order to fix a gauge, we should add to the action the gauge-fixed term 

1 



S. 



gf 



d'x Tr (A^V^)(V.A,) = -- / d^x (V^M;:, 



After some transformations the QCD action is 

+ J d^x Tr hgP^Mf^, A,] + |-[4, ij2 - 2B^,bA , 

The part of action quadratic over the quantum fields A" can be reduced 
to the form 



/ 



d x Tr 



(A^V.)(V.A^) + 8tgA,B^,A, = -^ I dV A"(G-^)"t^ 



b 



where 



(G- 



-l\ab 



-Vj5^^ - MQb^^ 



ab 



Finally, the QCD action reads 

So = -\jdxAl{G~XX = -\{\G-,lA. 
Si = 



(16) 



(17) 



9^ 



d^x Tr I igTf^^A^, A,] + ^[A^, A,f 



4 The Gluon Green function 

Let us find the gluon Green function corresponding to the operator 
The 8x8 matrix Ct for any vector n" has at least two zeroth eigenvalues. 
The choice of a vector n" requires special consideration and we plan to do 
it in further publications. In this paper, we choose the simplest case n^ = 
1, n"' = for a = 1,...,7 which guarantees the quark confinement (see 
[|, |, 0). In this case, it is convenient to introduce the notion 



e = diagd, 1,1, 0,0, 0,0,1 



, J-, J-, w, vj, vy, vy, -1-y, 



/i = diag(0,0,0, 1,1, 1,1,0), 
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Q 



ab 



e + h = I = diag(l, 1, 1, 1, 1, 1, 1, 1), 
^'''"' gB = A'V3, 



gB{-iri = ^n' 



ab 



n 



h. 



In these notion the operator ([l6|) can be written in the form 

{G-'),. = -dX'J + (-VJ + 2A')PM + (_vj - 2A2)P(;) (18) 
where the projection operators are introduced 



P,^? = e5, 



III/ 



liui 



p(±) = - 



6^,yh ± ihbi,„ 



(p(±)]T ^ p(±) 



Ufl 



Here the symbol " *^ " is connected with color indices. 
The equation on the Green function looks like 



(G-XuKix) 



'^l5,u 



AiVLh 



fii/ 



ab 



Kix) 



(19) 



-dX'J + (-VJ + 2A^)PW + (-VJ - 2A^)P(- 



Then, the Green function can be written in the form 
G^„y{x, y) = F^^{x -y) + D^^{x, y), 



Af,{x) = J^{x) 



F,,{x-y) = F{Z)P;^'J, 



F{Z) 



p 



6{z), 



(20) 



D^,ix,y) = e'^^^'y^D, 



D^,,{z) 



D^^\Z) 



-Vliz) ± 2A2 



{-) 
6{z) 



A222 



where z = x — y and Z ^ . 

The free part of the Green function is 



F{Z) 



/A 



:m'^'^-[Tj z 



(21) 
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The regular Green function D^^^ reads 



D^+\Z) 



1 


+ 2A2 V47ry 


/ du ■ ■ e 

J u+ I 
1 


/A\2 

\4:7tJ 


+ 2e^Ei( 2Z) 

Zj 


^ UttJ 2Z2 



-Zu 



for 



(22) 

Z ^ CXD. 



It is convenient to represent 

fl«(z) ^ (A) /,P< 



+) e-^'. 



(23) 



The operator — Vp(-2) — 2A^ in the function D^~\Z) contains the so called 
zeroth mode, i.e., there exists the eigen function $o for which the eigen value 
equals 



[-V^(^)-2A2]$o(^) = 0, $ 



ol2;J 



C^-e 



_A^ ,2 



In order to find the contribution of this zeroth mode to the Green function, 
we should solve the equation 

[-Vl{z) - 2X^]m,{z) = Jo ■ %{z) = Jo ■ e-^^^ 
This equation can be solved and the behaviour of this solution is 



"^niz) ~ e^""^ 



for 



oo. 



It means that this function does not belong to the space L^ and can not 
describe any physical state. Therefore, this function should be excluded 
from the set of physical states and the gluon Green function can not contain 
this contribution. 

Thus, the Green function D^~^ without the zeroth mode looks like 



D^-\Z) 



reg 



1 



-V2(^) 



2A2 



6iz) 
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where reg means that the zeroth mode should be excluded. For this aim let 
us use the resolvent 



1 

R(a = -^ 5(2)= /rfae"^'(")+(2A2-C)«^/ X 

-V2(z) - 2A2 + C ' ' / 

\ 2 



da e ^'^ 



2 ^A^c 



A^e 



^47rsinh(A^a;) 



e 4 



^^^^ coth(A2a) 



Uvrj 



du 



(coth(M) + i)2e-^^°th(«) _ 4g-z 



+ 



4A2 



-.-^ 



for C ^ 0. 

Then, the Green function D^^^ is 



D^-\Z) 



A \ ^ 7 (iw 



An J J u^ -I 
1 



47r/ 



[u + l)2e-^" - 4e-^ 



(24) 



2C-21n(2Z) 



■ e 



A 



A-kJ 



\ \n{2Zy 



for Z -^ 00. 



Here, C = 0.577215... is the Euler constant. We shall use the representation 



^'"'<^' = (s) /''p^'-"" 



(25) 



Finally, the gluon Green function reads 



G,,{x,y) = F{Z) ■ Pj^j + e^^^^'y^D,,{Z), 
D,.{Z) = D^i^z) . P(+) + D(-)(Z) . P/-) 



fiu • 



(26) 



4.1 Confinement of gluons. 

Confinement problem: four color degrees of freedom are confined, but the 
other four color degrees of freedom are not confined. Our assumption is 
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that the gluon vacuum has a domain structure R = IJ F^ with the scale of 

3 

distances 

J^hadron *%^ -^domain *%^ ^exp- 

Here Rhadron ~ ;!; is a scale of regions where the strong quark-gluon-hadron 
interactions take place, Rdomain is a size of domains Tj, R^xp is a macro- 
scopic scale where hadrons are registered. Each domain is characterized by a 
direction of the color vector n" and this direction is supposed to be random. 
Now we would like to give qualitative and semi-quantitative arguments 
that this structure of gluon vacuum can provide for the total gluon confine- 
ment if directions of the color vector rf" are random in each domain. Let 
us consider for simplicity the case of group SU{2). In this case we have the 
same problem. The 3x3 matrix n looks like 

^ab ^ _i^abc^c^ ^^c^c-^ ^ ^^ ^^^ b,C=l, 2, 3) 

and the eigen values are Aq = 0, A± = ±1. The qualitative behaviour of the 
gluon Green function is 

where Go is the free Green function. For n = (0, 0, 1) we have 

G^\z, n) = G^^z) = e-^, G^%z) = Go, 

i.e., we have confinement for directions 1 and 2 and no confinement for di- 
rection 3 in the color space. 

Now let us consider a gluon which moves in a direction through A^ do- 
mains. It is clear that the terms with e~^ provide the confinement in each 
domain but the terms with Go permit gluons to travel at large distances. 
Let a gluon moves in a direction through A^ domains. Its amplitude will be 
proportional to the factor 

Afi = nj(n/n7v)(n7vnAr_i) ■ ... ■ (n2ni)(ninj)nj. 

Here Uj is the color vector of the domain Tj. We can consider {rijUj^i) = 
cos{6j) where < ^j < vr is assumed to be a random angle. Then, one gets 



N ( ^/2 

^/il ^ IT I cos(^^)| = exp i iV f d9 \ncos{e) 



^-Nc 
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where c = 1.0887... . Thus, we see that a large number of domains can 
confine gluons. 



5 Two-gluon scalar colorless current 

The matrix element of the second order, which contains the scalar two-gluon 
current, looks like 



iw) 



2 JJdxdy T;,{x){-V;\x)V::{y)G%{x,y))T:^^{y) 



2 

{2t9y 



2/ ffdxdy T;,(x)fe^^(^'^)5^,,,/,(x-|/)e^^(^'^)l"Vj,(i/) 



where 



t;.(x) 



Tr ([i^(x),i,(x)r). 



and 



^S,./3(^) 



T,lix) 



Tr (i^(x)i,(x)r 



[V,{z)VMD^(^{^)r, 



(27) 



^taMlsi^) — '^uaMlsi^) + '^tl3,nai^) ~ ^"a,u/3(^) 



'^fi/3,ua\^>- 



To obtain the function ^^qj,^(-2) the following calculations must be ful- 
filed: 



yab/ 



-'^■'{^)'^a{y)G%{x,y) 



yab/ 



-V'!^{x)G%{x,y)Vr{y) 



-Vf{x)Vl''^iy)GT,{x,y) = -Vfix)Vl''^{y) [/'^^^'^D^,{x - y) 
- e*(^'^)V^(x - y)Vliy - x)D,p{x - y)]" 



Using the identity 



^aYg|b(x,j/)r'' ^ / lb(x,j/)^feg-lfe{:.,j/)\ 



one can obtain 

ab 



Tr {A^{x)A,{x)e) (et^(^'S')) = Tr {A,{x,y)A,{x,y)t'), 
(ei^'^^'y)YTT{fA^{y)Ap{y)) = Tr {fA^{y,x)Ap{y,x)) 
where the bilocal gluon fields have the form 

A^{x, y) = e-^^(^'^)i^(x)e^^(^'^\ A^{y, x) = e-5%.^)i^(y)ei^(?''^) 



So we can rewrite L2 as the follows 



L2 = 2g^ JJdxdy Tr(i^(x, y)A,{x, y)e)Sll,p{x - y)1i[fA^{y, x)Ap{y, 



In order to extract the scalar colorless two-gluon state, the following stan- 
dard steps should be done: 

1. Using the Fierz transformation for the t"-matrices 

we extract the four-gluon color singlet 



Tr [A^{x,y)A,{x,y)r)Sll,p{x-y)i:i [fA^{y,x)Ap{y,x) 
= Y^Tr (i^(a;, y)A^{x, y)Aa{y, x)Af^{y, a;))5^^_^^(x -y) + .... 

2. Next we extract the two-gluon color singlet and octet: 
Tr (^Af,{x,y)A:,{x,y)Aa{y,x)A^{y,x) 
= (Au{x, y)Aa{y, x)) . . (A^{y, x)A^{x, y) 

= o^r [K{x,y)A^{y,x)^ -Tr {^Afi{y,x)A^{x,y) 



+2Tr [eA,{x,y)A^{y,x))-Ti [eAp{y,x)A^{x,y 
3. The two-gluon color singlet can be represented 

J>.a{x,y) = ^J{x,y) + jl^Xx,y), 
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where J'{x,y) and jj}^{x^y) describes spins zero and two respectively. For 
spin zero we have 

J{x,y) = 2Tr (i.(x,t/)i,(i/,x))=2Tr (i.(x)e^^(^'^)i,(y)e*^(^'^); 

= Al{x)(e'^^-'y^Y' Al{y) = (A.(x)e*^(^'^M,(|/)) . 

4. The effective "potential" connecting two scalar currents is 
1 



V(^) 



Finally, we get 



" oaa 

7^ fiV,U^ 



(^) 






(^) 



(29) 



Tr 



9' 



V,{z)Vl{z)D,,{z) - V,{z)Vl{z)D, 



L2 = — // dXidX2 J{xi,X2)V{xi - X2)J{x2,Xi) 

5.1 The gluon-gluon "potential". 

The gluon-gluon "potential" V{z) can be calculated 

V(^) = Tv[v,{z)Vl{z)D,,{z)-V,{z)Vl{z)DU^). 

= Tr iv,(z)pWvJ(^)-V,(^)V;[(^)pW]DW(Z) 

V,(.)V^(.)P(;)Ld(-)(Z). 



(30) 



Tr 



v,(^)p(;)vj(.) 



One can get 



VJ:)P'.i>^J(z)D'^>(Z) - V„(2)VT(2)P»>B'±>(Z) 



"i 



-V^(z)±2A^' 



T 



T^-6Z±2Z— 
dZ 



D^^\Z)\. 



The first term gives 5{z) and this part of our potential describes the local 
interaction which leads to the standard renormalization and contributions to 
the highest perturbation orders. Thus, the potential that is responsible for 
the glueball bound state reads 



W{z) 



4A2 



3Z 



dZ 



D^+\Z) + 



2-3Z-Z 



d 

Iz 



2^ 



2 00 



dWa e 



~Za 



D(-)(Z)| 

(31) 
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It is usefull to introduce the representation 

A 



W{Z) = — jdUae-^\ (32) 



6 Nonlocal Currents and Bosonization 

Now let us go to L2 in (^) and consider the scalar current J{x,y), which in 
the center of mass system of two gluons (xi = a; + |, X2 = x — ^j becomes 

y 

X 

2 



\(x)ei^("M,(x; 



where 



V/. {x) =V/. (x)- V,, (x), 

V/. (x) =d,i -iQbf,„Xy, V/i {x) =8,1 +iVLh^yXy. 



Then 



L2 = ^JJ dxdy J{x,y)Wiy)J{x,-y) 



dx dy,dy2^JW{yi)J{x, yi)5\yr - y2)\/v%^J(a;, -y^) 



= l- 

864 

Let us introduce an orthonormal system of functions {Uq{ij)\ 

E f^Q(2/i)f^Q(2/2) = 5(1/1 - 1/2), \dy \JQ{y)\JQ,{y) = 6qq,. (33) 

where Q are quantum numbers. We obtain 

2 
L2 = ^Jljdx Jq{x) . J^{x) (34) 

w 

where 

Jq{x) = (^A,{x)VQ(y {x))A,{x) 

Vq{V (x)) = jdy ^W(^)UQ{y)e^^^\ 
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We are interested in the scalar glueball which is described by the current j7o, 
i.e., with the quantum numbers Q = 0. The part of L2 with these quantum 
numbers looks like 

L, = ^-^JdxJ,\x), gf = ^ (35) 

Jo{x) = (^A^{x)VoiV ix))A^{x] 

The scalar gluon current consists of the confined and non-confined gluon 
modes 

J = j(f) + j(-) = (4^VoA(/)) + (4^Vo4^)) 

It should be stressed that this division is SU{3) invariant because the eigen 
numbers of the matrix n are invariants of SU{3) group transformations. 
Therefore, 

(A.G-^lA,) = iAi^G-^lA^:\x)) + {AfG-lA\J\x)). 

Our main assumption is that only confined modes lead to a bound state, 
but local modes lead to standard renormalizations and contributions to next 
perturbation orders. It means that the part of the partition function, which 
is responsible for the glueball formation, reads 



L)V9e~5(w)-5Trln(l-x/29i'pV'oD)^ 

Expanding the logarithm, we extract the quadratic part over '^{x) in the 
form 

where 

H~^{xi - X2) = 6{xi - X2) - gJTl{xi - X2), 
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U{xi-X2) = Tr {VoiV.c^)D{x,,X2)Vo{Vx,)D{x2,xi)]. 
The polarization operator Il{z) can be written as follows {z = xi — X2) 

U{z) = jj dy^dy2 v/W(i/i)f/o(yi)n(^||/i,t/2)v'W(i/2)f/o(y2) (36) 

where 

^{z\yi,y2) 

= Tr je^^("^)e*^("^^"^)D«/3(a:i - X2)e^^("2)e*^("^^"^)Dj„(a;2 - xi) 

Tr je^^«^-^^)^^)Z}„. f^ - y-^^\ Dj^ ( -z ^^ + ^' 



^a/3 



2Tr 






with 



A2 



^(±) 



2± 



yi+1/2 



6.1 Polarization operator 

We do not know the eigen function Uo{y); so we shall use the variation 
method with the test function 

A^ A2 

Ztt 4 

where 7 is a variation parameter. Our point of view is that this function 
should be a quite good approximation for the ground state of the Bethe- 
Salpeter equation. 

The Fourier transfrmation of the polarization operator looks like 



n(p2, 7) = / dz e'P'U{z 






R{p\i) 



(37) 



R{p\ 7) = 7' jj dUa.dUa, E jj dV^^UV^J I{a,, 02, c, + 7, C2 + 7; p')- 
1 i=± 1 
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/(ai,a2,Ci,C2;p^) = Ti 2h dz jj dyidy2 

^pz-a^^-^~a2'-^~c^4(z~!ilp^f~c,4(z+y-lp^)\in((y^-y2)bz) 



6 O „2 



47r\ 8 p2 fc(a^,a^,ei,e2) 
I . . g A^ if(a]^,a2,ci,C2) 

AV [ir(ai,a2,Ci,C2)]2 
A;(ai, 02, ci, C2) = 0102 + -(oi + a2)(ci + C2), 
K(ai, 02, ci, C2) = (oi + a2)(l + C1C2) + (ci + C2)(l + 0102). 

Finally, we have 

9IW.1) (38) 

2 00 00 { P^ fc(<iiiC2,ci+7,C2+7) \ 



TT 27 77 '•^ ,ti77 ^1 ^^ [ir(ai,a2,Ci + 7,C2 + 7)] 

1 ''I 



2 



where a^ = ^. Next step is to calculate 



2^( M^ \ as ^(M^\ 
and the mass of scalar glueball is defined by the equation 

2\ 



^-(^1- (-) 



The results of calculations are given in the Table. 
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M 
A 


as 


M 
(A = 446 Mev) 


3.0 


1.5120 


1337 


3.1 


1.0609 


1381 


3.2 


0.7369 


1426 


3.3 


0.5067 


1470 


3.4 


0.3450 


1515 


3.5 


0.2325 


1560 


3.6 


0.1551 


1604 


3.7 


0.1024 


1649 


3.8 


0.0669 


1694 


3.9 


0.0432 


1738 


4.0 


0.0277 


1782 



The gluon condensate 

\ TDO, DO. 

ill 

1 






rya rya 
^^iv^fiv 



nVB\,b^, 



AB' 



2 d2 _ ^ y . ^q2 



A4 = Zg^B'^ = — 

3 3 47r2 



TT^ a.c 



-^■^■{G;^Gl^) = -{Oa) 



TX 



(Og) = M2xG {Gev)\ 
A = U6xf Mev, 

where xg ~ 1 (see, for example, 0]). The values of as(A^) are given in [|ll 



6.2 Conclusion 

In the conclusion we would like to list the main points of our approach. 

• Analytic confinement plus weak coupling lead to the Regge spectrum 
of bound states. 
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• Gluon self-dual vacuum field B^{x) leads to full quark confinement and 
partial gluon confinement. 

• The field B^{x) has directions in color and configuration spaces; so the 
domain structure of the whole space R^ with appropriate averaging 
over these directions is needed to provide full gluon confinement and 
relativistic covariance. 

• Gluon propagator in the field B^{x) has the zeroth modes, which are not 
the Goldstone degrees of freedom. We exclude them, because they lead 
to a nonphysical growth of the propagator in the x-space for x ^ oo . 

• Numerical results are quite reasonable and with the results of [^ ^ 
give hope that our main assumption about B^j_[x) contains "a grain of 
truth" . 

6.3 Appendix 

We use the following method of approximations. Let a positive function 



f{t)=Jdaue-'^ (40) 



be given, where dau is a positive measure. 
Our aim is to get the approximation 



f{t)^Ut) = Y.A,e-^^' (41) 



on the interval 0<ti<t<t/<oo. 

Our procedure is the following. Let us choose n points 



and calculate 



{tl, ^2, ..., t„}, tj G [ti,tf\ 
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Then, we solve the system of equations 



/(t,) = ^A, e-«^*% (^ = l,...,n) 



so that we get the numbers 



{Ai,...,An} 



—a-jt 



and we have the approximation 

oo 

f{t)= /da„e-*"^/„(t) = X:^, e-" 
1 J=i 

The accuracy of this approximation can be evaluated by 

llf - f II / 

11/ /all M.M_ / ^^1^(^)1 



(42) 



The accuracy depends on a particular choice of the points {ti, ...,t„}. 
In our case, we have 













D(-\t) 


V'w(t) 


D(+)(t) 


t, 


.04, .4, 


1, 2.2, 4, 7, 10 


.01, .05, .2, .4, 1.2 


.01, .05, .1, .3, .5, .9, 3 


j 


^. 


Cj 


^. 


ttj 


^. 


ttj 


1 


7.009 


3.341 


179.4 


106.3 


175.2 


96.62 


2 


0.485 


2.227 


29.42 


24.04 


36.16 


19.76 


3 


-8.234 


0.959 


6.334 


7.607 


152.7 


10.94 


4 


-1.532 


0.002 


0.281 


4.592 


27.77 


6.554 


5 


55.497 


0.144 


1.031 


2.365 


-171.8 


10.46 


6 


-257.1 


0.197 






3.971 


1.836 


7 


208.3 


0.215 






-2.783 


0.823 


5 


0.041 


0.027 


0.035 
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